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^ ■ Abstract 



O \ We study the reconstruction of the source function in space-time directly 

^3 



from the measured HBT correlation function using the Maximum Entropy 



^ \ Principle. We find that the problem is ill-defined without at least one addi- 

H : 

>y , tional theoretical constraint as input. Using the requirement of a finite source 

lifetime for the latter we find a new Gaussian parametrization of the source 
function directly in terms of the measured HBT radius parameters and its 
lifetime, where the latter is a free parameter which is not directly measurable 
by HBT. We discuss the implications of our results for the remaining freedom 
in building source models consistent with a given set of measured HBT radius 
parameters. 
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Relativistic heavy ion collisions are used to create small and shortlived regions of hot 
and dense hadronic matter under conditions resembling those of the Early Universe near the 
confinement phase transition. These little "fireballs" expand rapidly into the vacuum, finally 
decaying into a multitude of particles whose momentum distributions can be measured exper- 
imentally. The only known practical way to extract direct information about the space-time 
distribution of the source from such momentum measurements is via Hanbury-Brown/Twiss 
(HBT) intensity interferometry which exploits the quantum statistical correlations in the 
2-particle coincidence spectra. The latter are reflected in the well-known relationship be- 
tween the momentum correlation function C {q, K) of two identical particles and the source 
distribution function S{x, K) 



Here, K = (pi + p2)/2, = {E^ + E^)/2 ^ Ek = Vm^ + K^, q = pi - p2, g° = - E2. 
Our interest focusses on the source function S{x, K) which describes the probability that a 
particle with 4-momentum p is emitted from the source at spacetime point x. It determines 
the effective spatial and temporal size of the region emitting particles with momentum 
K. The two-particle HBT correlation function C{q,K) is usually well parametrized by a 
Gaussian function in g^=(go5 Is, Qo, Qi) H] where go, li, lo, Is denote the components in time 
direction, parallel to the beam ("longitudinal" or 2;-direction), parallel to the transverse 
components of K ("out" or x-direction), and in the remaining third cartesian direction 
("side" or y-direction), respectively. For azimuthally symmetric systems it can be written 
in the "standard form" as 

C(q, K) ^ 1 ± exp [-qlRl{K) - qlRl{K) - qfR^^) - 2qiq,Rl{^)\ , (2) 

or in the so-called "Yano-Koonin-Podgoretskii (YKP)" form as 

C(q, K) ^ 1 ± exp [-gli?l(K) - {q\ - ql)R\{Y.) - (g-t/(K))2(i?^(K) + i?J(K))] , (3) 

where 



[/(K)=7(K)(l, 0, 0,t;(K)), ^ = \l^fY^, (4) 

and qj^=\Jql + ql- Please note that since the momenta pi^2 are on-shell, the four components 
are not independent: they satisfy K ■ q = Q oi 

go = /9-q with {3 = ¥./ ^ K/ Ek ■ (5) 

(3 is the velocity of the particle pair. The mass-shell constraint (|]) is the reason that, 
for azimuthally symmetric systems and in the Gaussian approximation P|, the correlation 
function is completely determined by only four i^-dependent functions, the HBT radius 
parameters -Rf(K) ( i = o,s,l,ol) in the standard parametrizations, or -Rf(K) ( i =± 
,11,0) together with v{K) in the YKP parametrization. As discussed in Refs. P,||J^ these 
measurable parameters are linear combinations of the 7 independent second order space-time 
variances (rms widths) of the source function S{x,K). Thus the HBT correlator provides 
insufficient information to reconstruct the spatial and temporal structure of the source, even 
in the Gaussian approximation. A separation of the spatial and temporal structure of the 
source is thus only possible with the help of specific model functions for the source. So far 
the range of possible or even reasonable models is poorly understood, and its exploration 
remains an important problem in the context of HBT interferometry for the practitioner. 

In this note we want to try another approach, i.e. to construct a "most likely" or 
"least biased" class of source functions starting directly from the (insufficient) experimental 
information in terms of the HBT radius parameters, adding a minimal amount of theoretical 
prejudice in the form of additional constraints. It is well known that the Maximum Entropy 
Principle (MEP) [jlO[ is a powerful method in solving this kind of problem. We will study how 
the additional theoretical constraints affect the shape of this "most likely" source function, 
and to what extent such an approach can be used with advantage in the interpretation of 
HBT correlation data. 

As is well known, for a given normalized probability distribution p{x), the information 



entropy is defined as |]10 
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p{x) Inp(x) d^x . (6) 



In our situation the role of p{x) is taken by the source function S{x, K), where K is to be 
treated as an additional parameter. We are interested in the "most likely" spacetime distri- 
bution Sk{x) = S{x,K) for each value of K. Correspondingly, Eq. (||) must be generahzed 
as 

aK = - J S{x,K)\nS{x,K)d^x. (7) 
Let us assume that we have experimentally measurable constraints: 

{g,{x)) = G,{K) , z = l,...,N, (8) 

where the expectation value is to be taken with the source function: 

_ Jd^xgiix)S{x,K) 

In the case of HBT interferometry the functions gi{x) are certain linear combinations of 
bilinears in the components of the coordinate x^. The "most likely" distribution function 
S{x, K), which has least prejudice or bias, is the one which maximizes the entropy ax under 
the constraints (^). Using the method of Lagrangian multipliers one obtains 



S{x, K) = exp 



N 



(10) 



-\iK)-Y.KiK)g,ix 

i=l 

where A(K) ensures the normalization of S{x,K): 

j d^xS{x,K) = l =^ = J d'^xe-^^^'^'^^'K (11) 

(Usually the source function is normalized to the single particle spectrum, / d'^x S{x, K) = 
E dN / d^K, but the resulting modifications are trivial.) Eq. (|TT|) expresses A(K) as a function 
of the other Lagrangian multipliers Ai(K). These in turn are determined by inserting the 
actual values of the constraints into 

- = {g^{x)) = G,{K) , z = l,...,N, (12) 



which follows directly from Eqs. (|§), (J^) and ([lOD. 

It is instructive to first consider a simple though unrealistic case. If the HBT correlator 
would simply measure the relevant geometrical "radii" of the system, i.e. if 



GsiK) 



if) 



(13) 



where 



(/i = 0,l,2,3. 



(14) 



is the deviation from the (K-dependent) point of maximum emission, a;''(K) = {x'^), such 
that Eqs. (|T3|) define the variances of the source function S{x, K) in the various space-time 
directions, then the "most likely" source distribution (|iy) simply becomes 



S{x,K) = \{ 



1 



exp 



(15) 



=0 pT^RKK) 

This describes a static Gaussian source with i^'-dependent radii Ri{K), i = 1,2, 3, and finite 
/T-dependent lifetime Rq{K). (Remember that K is the momentum of the observed particle 
pair.) 

More generally, if the experimentally measurable constraints Eq. (^ are generated by 
bilinear functions of x. 



9ii^) -^/i 9i -^u 1 ^ I5 • • • 5 1 



(16) 



the resulting source distribution will be a general Gaussian, 



S{x, K) = e-^(^) exp 



N 



^ Ai(K) x^ g 



X]_j 



i=l 



e-^(^) exp 



x^B^''{K)x, 



(17) 



which agrees in its form with the generic Gaussian saddle point approximation |^,[n]] of 
S{x, K) around x(K). Here, B^''' = 2Y.f=i Kg^ and e'^ = VdetB /{4n'^). B is a symmetric 



4x4 matrix and has usually 10 independent parameters. For an azimuthally symmetric 
source, only seven of these are non-zero 0. 

Now let us turn to the situation as it is actually encountered in nature. If the measured 
correlator C(q, K) is fitted only to the standard form (0) we know that the extracted HBT 
radii correspond to the following constraints on the source [§,0: 

Ri = m , 

Rl = {{x-(3otf), 

Rl = {{x-(3oi){z-f3ii)). (18) 

In this case the normalization constant is found to be e'^ = {27r)^^'^ Rs^J R^Rf — {R^iY I di 
and thus diverges. Furthermore, the corresponding matrix 

/ 7-)2 n r-)2 Ty2 o o \ 



1 

A 





0, 


-Rill 


RlA 


- RWl 


0, 




0, 







-Rlh 


0, 


Rli 


Rlif^i 


-R% 



(19) 



Ri,(3i - RfP^, 0, RiiP^ - RiPi, RfPi + Ri(3t - 2i?^,/3^A ) 



with 



2 \2 



A = RtRt - (Ri,) 



(20) 



has a zero determinant. Modulo the diverging normalization this implies that the t, x and 



z variances of the corresponding source function ([T^) diverge. All this implies that the 
measured HBT radii (|18]) are not sufficient to constrain the space-time structure of the 
source in all 4 space-time directions. The reason for this is that the observed particles are 
on the mass-shell and that therefore, as mentioned before, the time-component is not an 
independent variable which can be used to explore freely the necessary fourth dimension. 
Thus, based on purely experimental information from the standard fit alone, it is unable 



to obtain a well-defined four- dimensional space-time distribution of the source. We have 
checked that one arrives at the same conclusion starting from the YKP fit instead of the 
standard fit (Q). This is to be expected because the two fits are mathematically equivalent 
and correspond only a different way of eliminating the redundant component of using 
Eq. (^). The similar problems encountered with the Maximum Entropy Method in both 
cases reflect this equivalence. 

It is clear from this discussion that the Maximum Entropy Methods tends to construct 
from the data a source with an infinite lifetime. This is certainly not reasonable for the 
sources created in heavy-ion collisions, since they disperse very rapidly into the surrounding 
vacuum. But even in other situations, where it may not be a priori clear that the source 
lifetime is finite, one can easily check with the available experimental information whether 
the source is stationary or not. For a stationary source one sees from Eq. (|l]) that the 
correlator can only be different from unity if g'' = /? ■ q = 0. This means that for pairs at 
y = = there cannot be any correlations in the "out" direction; in other words Rq 
must be infinite. A finite value for Ro is thus not compatible with stationarity of the source. 
In this case, in particular for heavy-ion collisions, it is thus natural to require, in addition 
to the experimental constraints (|^), the finiteness of the source lifetime: 

(?) = T\K). (21) 

Here T{K) is an unknown, but finite A'-dependent parameter. We will now show that this 
leads to a well-defined source function. It will be parametrized by the chosen function T{K), 
in addition to the measured HBT radii. 

Adding Eq. (^) to the constraints (§) and going again through the procedure of maxi- 
mizing the entropy ax we now find 



exp 

S{x,K) = 



' 2Rl 2A 2T2 



[22] 



Ati^TRsVA 

with A from Eq. (PP]). Note that this source function depends, via the definition ( ]T^ ) of 
x^^, on four arbitrary, K-dependent parameters x'^(K) which determine the position of its 



saddle point (a;'^)(K). The latter cannot be determined by measuring particle momenta 
since an arbitrary (even i^'-dependent) translation of the source does not influence the latter. 
Thus any source function which is able to reproduce the measured momentum spectra and 
correlations will have this remaining ambiguity. 

It is interesting to study the general physical properties of the source function in 
more detail. The corresponding matrix B is now regular and given by 

/ T-.9 ^ 7-.9 7-.9 o 7-.9 ^ ^ 



1 

A 





0, 


-R-lh 


R-lA 




0, 


A/Rl 


0, 







-Rlh 


0, 






-RlPi 



[ RlA - Rf(3^, 0, - Rl/3i, Rf/3l + Rl/3f - 2Rli(3^/3i + A/T^ ) 

(23) 

It posesses the following important properties: flrst, it has the correct symmetries in the 
sense that in the limit — > (where the x and y coordinates are indistinguishable and 
therefore Ro = Rs = R±, Roi = 0, and A = R\Rf |^), the terms -B13 and Bio vanish: 



lim (Bf,^) 



( 1 

0, 
0, 
0, 



0, 0, 



0, 



0, 
1 






Pi 



n -A ^' -^ J- 



1 

0, 
0, 
0, 



0, 
1 

0, 
0, 



0, 

0, 
1 








(24) 



Secondly, for nonzero values of the nondiagonal elements of B generating x-z and x-t 
correlations do not vanish. In contrast to the models studied in Ref. there is actually 
no reason for them to be even small. The nonvanishing components of the inverse matrix, 
h = (B^^), are given by 



x^)=bn = Ri + (3lT 



2 rTi2 



if) = b 



22 



{~Z') = &33 = Rf + f^fT' 

{P) = boo = , 
{~zi)=ho = PiT\ 
{xi) = ho = PlT^ , 



(25) 



Obviously, all second order variances are now finite. Also, one easily checks that in the limit 
^ the last two components, 610 and 613, vanish as required by azimuthal symmetry. 



The source (|2^) also gives the correct YKP fit parameters. Calculating the latter from 



with {C, = X + iy) 



Rn 



A 
B 
C 



= if) , 
A + B 
2C 

B-vC, 
A-vC, 



2C 



A + B 



t 



(26) 



(27) 
(28) 
(29) 



evaluated with the variances (pSf ) calculated from the source (p2D, one easily verifies that 
they satisfy the relations with the standard parameters ( [T8| ) given in Eqs. (21,22) of Ref. [^. 
However, since the experimentally unmeasurable parameter T from the additional theoretical 
constraint (^T]) is specified in a fixed frame, while the YKP parameter Ro measures the source 
lifetime in the comoving Yano-Koonin {v = 0) frame, there is no simple relation between Ro 
and T. For the same reason, the T-dependence of the space-time variances ( ^51) will look 
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less simple in any other frame which moves with a fixed, but non-zero longitudinal velocity 
relative to the frame in which (|2T| ) is specified. 

It is interesting to compare the matrix ( pSf ) with the one extracted in Ref. [^] from the 
constraints (|18|) under the additional assumption that the x-z and x-t correlations B^z and 
B^t are small and can be neglected: 



/ 1 



(5, 



0, 



n 



V 



0, 
0, 



0, 
0, 



0, 

0, 

Rl-Rl 






{R;i-Ri)Rf-(R^J2 



\ 



(Rfimw^iKi) 

{Rl-Rl)Pi-RliPi_ RfPl+iRl-Rs)Pf-'2Rli(^±Pi 



(30) 



This has already the block-diagonal shape (^41) , even for no n- zero (3±. It is a regular matrix, 
with non-zero determinant: 



detS 



(31) 



Rt mi - Rl)Rl - R%] ■ 

This means that the above additional assumptions were also sufficient to guarantee a well- 
defined and finite source function in the Gaussian approximation, using only the measurable 
constraints (p^Sf). These assumptions can be justified for sources without collective transverse 
expansion, where they are satisfied exactly |]5|jr3[|. For non-zero transverse fiow, however, 
they are usually violated P,p!3|. Our solution (p3|) does not make such assumptions, at the 
expense of introducing an additional unknown parameter function T{K). It is interesting 
to note that in the limit — the matrix (BOl) takes a form which differs from 



il^o ^^^^^ 





0, 


0, 





0, 


1 


0, 





0, 


0, 


(P) 


(t2>(52)-(t5)2 


0, 


0, 


m 


(52) 




{P)(z'2)-(izy 



\ 



(32) 



(Of course, (i^) = (y^) = R]_ in this limit.) The corresponding source (0) thus differs from 
the P± ^ limit of (^). Still, both sources reproduce the same constraints (|T^. This 
means that along the particular space-time directions in which the HBT correlator is able to 



10 



probe the source S{x, K) via the restricted (see Eq. (^) Fourier transform (|T]), both sources 
have the same Gaussian curvature. 

In this note we have taken the basic philosophy of HBT interferometry, namely to ob- 
tain information on the space-time structure of the source from momentum measurements 
only, to the limit. With the help of the Maximum Entropy Principle we derived the source 
S{x, K) directly from the experimental constraints provided by the measured HBT correla- 
tor, supplemented by one minimal but necessary additional theoretical constraint limiting 
the source lifetime. This additional constraint is required to render the problem well-defined, 
because the kinematic restrictions imposed by the mass-shell constraint @ generally pro- 
hibit the inversion of the relation (|l|) between S{x, K) and C(q, K). Our resulting source is 
a Gaussian in the space-time coordinates, parametrized by the K-dependent HBT radii (0) 
and the lifetime parameter T(K), as well as by the (unmeasurable) K-dependent position of 
its saddle point (x^). It differs from the Gaussian approximation given in Ref. which was 
derived using different additional theoretical constraints, namely the vanishing of the x-t 
and x-z correlations (xt), {xz) independent of the pair momentum. Both sources reproduce, 
however, the same (measured) HBT correlation function. Our result thus provides not only 
a possibly useful new source parametrization, but more importantly an explicit and inter- 
esting example illustrating the remaining freedom in the space-time structure of the source 
left open by even the most accurate HBT measurements. 
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